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ABSTRACT 

I derive conservative upper bounds on the supersymmetry breaking parameter m 1( / 2 
as a function of the strong coupling in the Standard Supersymmetric Model (SSM) 
using gauge coupling unification. I find that over more than 99% of the parameter 
space, «3 > 0.120 implies that mi/ 2 is below lOTeV and 0:3 > 0.129 implies that 
mi/2 is below 1 TeV. I express the variation of these bounds over the SSM parameter 
space with a numerical coefficient, c. I also find that in the SSM, a reasonable value of 
50 GeV < mi/2 < 1 TeV requires a 3 > 0.119 over the whole parameter space. These 
bounds are particularly sensitive to the value of sin 2 9w = 0.2317 ± 0.0005 used in 
the calculation. In more realistic models, heavy thresholds and gravitational effects 
will modify this result. Although these effects are theoretically calculable in specific 
models, more realistic models contain many unknown parameters in practice. I illus- 
trate this point with minimal supersymmetric 577(5) where the combined constraints 
of gauge coupling unification and proton decay require a 3 > 0.119 for mu2 < 1 TeV 
and the upper bound on the supersymmetry breaking scale is greatly relaxed. 
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Gauge coupling unification JIJ applied to precision LEP measurements has pro- 
vided strong motivation for supersymmetry, || and inspired numerous attempts to 
extract the supersymmetry breaking scale and constrain the parameter space of spe- 
cific models [Q. However, definite conclusions only result from very specific models: 
threshold effects and gravitational corrections make model-independent statements 
difficult M. In this paper, I find interesting constraints from gauge coupling unifi- 
cation within the specific framework of the Standard Supersymmetric Model (SSM) 
as defined in reference || and the minimal supersymmetric SU(5) model. Although 
these models are an excellent beginning, they can hardly be considered ultimate theo- 
ries as neither include gravity and the fine-tuning problem in minimal supersymmetric 
SU(5) requires some modification of the GUT structure. However, there is hope that 
similar constraints could be derived for specific realistic string models where the ad- 
ditional threshold and gravitational effects are in principle calculable. 

Direct searches for supersymmetric particles have continually increased the 
lower bounds on their masses. However, the only upper bounds on supersymmetric 
masses come from naturalness arguments |J or cosmological constraints on the LSP 
relic density j7|. Although compelling, the naturalness bounds are not rigorous, and 
whether the upper bound on supersymmetric masses is 1 TeV, 10 TeV or even 100 TeV 
is not clear and somewhat a matter of taste. The cosmological bounds can be evaded, 
for example by breaking R-parity. It would be extremely useful to have some other 
method of bounding the supersymmetric masses. In this paper, I focus on a first step 
in this direction by deriving an upper bound on the soft supersymmetry breaking 
parameter mi/2 as a function of the strong coupling in the Standard Supersymmetric 
Model (SSM). Our approach attempts a general and analytic analysis to reveal the 
underlying physics, and is meant to complement the many numerical searches of the 
SSM parameter space in the literature ||. Unless otherwise indicated, our notation 
assumes all gauge couplings are renormalized at mz in the MS scheme. 

Ignoring gravitational effects, gauge coupling unification in the SSM gives a 
simple prediction for the soft supersymmetry breaking parameter m^. Restricting 
attention to values of mu2 for which all the supersymmetric thresholds are above mz, 
this prediction is ||: 

ln(^) = -X + ^ + 71n 

m z a 3 

where the gaugino masses are given by 

oiAni , i 

m l = c l m i/2 = mi/2 , (2) 

«g is the unified coupling at the unification scale, and f(y,w) includes the threshold 
effects of squarks and sleptons under simple assumptions for the form of the stop 
mass matrix || 
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f(v,w) = y ln ( v c 9 + y) - 4 ln (\M + y) + 2 ln (V c ^ + y) 
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where y = mj/m^ and w = m 2 /rriy 2 with m 2 representing the off-diagonal elements 
in the stop mass matrix. The quantity X is given by 



X 
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.2 + 5 s (gauge) + 5 s (Yukawa) — sin 2 6 
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where /i is the dimensionful higgs mixing term in the SSM superpotential and mis- 
represents the mass of the charged, pseudoscalar and heavy scalar higgs which are 
nearly degenerate in the SSM. The gauge numeric corrections to two loop accuracy 
are || 

5 s (gauge) = 0.00127 + 0.01480a 3 , (5) 



and numerical calculations give —0.0004 < 5 s (Yukawa) < QT0 1 . To derive the most 
conservative upper bound on mi/2, I will want to maximize the right side of Equation 
(1), and therefore minimize X. Note, by definition, there are no GUT or intermediate 
thresholds in the SSM §. 

The calculation is most sensitive to sin 2 6^, for which the central value has 
come down in the last year to JTT 



sin 2 9 W = 0.2317 ±0.0005 . (6) 
Bounds on other parameters I will use are: 

a em = 131 GeV < m t < 190 GeV . (7) 

127.9 ±.0.1 v ' 

Although the most conservative choice for \x and m# would be to take them 
below mzi values of these parameters in the SSM turn out to be nearly proportional 
to m\i<i because of the correlations introduced by radiative electroweak symmetry 
breaking. In fact, /i and tuh are determined as a function of the five parameters of 
the SSM: m t , tan/?, my2, £o, £a- The explicit form of the tree level expressions for \i 
and nriH show that they both approximately scale with m^: 

m l,H = a ^H m l/2 + b v,H m w ■ ( 8 ) 

I use a parameter, c, to encode the dependence on \i and tuh in a simple way by the 
equation 

4 m z ' K m z J 4 K m z J y ' 

Note that throughout this paper, the logs associated with a particular threshold 
correction are understood to be zero if the threshold is below mz- A Monte Carlo 
search of parameter space reveals that c > 0.5 for more than 99% of the points 
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considered. Our search uses mi/2 = 1 TeV and randomly selects 10, 000 points for 
each sign of fi over the four dimensional parameter space defined by 130 GeV < m t < 
190 GeV, 1 < tan/3 < 50, -10 < £ < 10, and -2|f | < U < 2|£ |- The parameters 
m t,£,o,£,A are searched on a linear scale while the parameter tan/3 is searched on a 
logarithmic scale. 

To reliably extract mi/ 2 , the dependence of the gaugino masses on mi/ 2 (the 



EGM effect) must be carefully taken into account |12| . To derive the most conservative 
upper bound on m 1 / 2 , I must use an upper bound on a 3 (rrig), and and a lower bound 
on a 2 (fHii). The simplest way to extract the gauge couplings renormalized at the 
gaugino mass is to iterate the expression 

oi ■ 

aAmi) = — -, — r (10) 

lK l) l- a Aln( ai(m >^ ) 1 J 

along with the expression for the gaugino mass, Equation (2), to a solution. Since 
these values for the gauge couplings increase with bi, I take the maximum value of 63 
below the gluino threshold and the minimum value of 62 below the wino threshold to 
derive the most conservative bounds. Since in the SSM, the squarks cannot be much 
lighter than the gluino, I use 63 = —7. Below the wino threshold, I use the minimum 
b 2 = —19/6. Numerical calculations show that an upper bound on «g gives the most 
conservative bound for m^. I obtain this upper bound on a G by extrapolating the 
hypercharge coupling a y {mz) to the scale 10 17 GeV using the maximum b y = 33/5. 
This results in «g < 0.0454 after numerically correcting for two-loop effects. 

Finding an upper bound on f(y, w) depends on the physically acceptable range 
of the variables y and w. Although the minimum value of f(y,w) = —0.025 has 



already been determined |L3|, the maximum value is not as clear cut. For fixed 
w, f(y,w) approaches zero as y becomes very large. Maximization with respect 
to values of w yielding positive m~ gives a maximum of f(y,w) = 35/n( v /y)/48 at 
w = Cg + y as y becomes very large. However, a more physically reasonable choice 
is \w\ < dm t /mi/2 where d is a positive numeric coefficient of order one. Using this 
form for w, I have numerically scanned over values of d < 10 to determine that, for 
mi/2 > 1 TeV, f(y, w) < 0.5, which I will use in the calculation. It is the fact that the 
SU (3), SU (2), U (1) beta functions are all equal for an entire generation coupled with 
the assumption of universal soft supersymmetry breaking which severely restricts the 
range of f(y,w) and makes the results of gauge coupling unification approximately 
independent of the universal scalar mass itlq. 

Putting all this together gives equations which can be iterated to find an upper 
bound on m^j as a function of the strong coupling. The choice of parameters I make 
to minimize X and thereby obtain a conservative bound are: 

a em = -^— m t = 131GeV f(y,w) = 0.5 a G = 0.0454 . (11) 
LZ 1.0 

The bound using c = 0.5 is plotted as solid lines in Figure 1 for the central and 1 — a 
values of sin 2 6V- If «3 > 0.120, m 1 / 2 is constrained to be less than 10 TeV and if 
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ct 3 > 0.129, mi/2 is constrained to be less than 1 TeV. These solid lines would move 
downward as c increases giving tighter bounds on m 1( / 2 and would move upward as c 
decreases giving looser bounds on mi/2- To quantify this effect, in Figure 2 I indicate 
the value of a^, as a function of c, above which 777.1/2 is bounded by 1 TeV, 10 TeV 
with solid, dashed lines. In our Monte Carlo, among the 12,231 out of 20,000 points 
which have perturbative Yukawas and a stable electoweak breaking minimum, the 
values of c are distributed as follows: 62.9% give c > 2, 29.9% give 2 > c > 1, 6.4% 
give 1 > c > 0.5, 0.7% give 0.5 > c > 0.2, and 0.1% give c < 0.2. The points with 
c < 0.5 all have /i/mi/2 < 0.4 but rarely have < 1- This observation and 

the form of Equation (9) indicate that the regions with small c are the regions where 
/V m i/2 is fine-tuned to be small. 

Previous work pursued an alternative to finding an upper bound on the su- 
persymmetry breaking scale and derived a lower bound on a 3 as a function of m^j 



13] . For values of 777.1/2 giving wino masses above mz this derivation followed from 



Equation (1). For values of 777.1/2 giving wino masses below mz, the relation becomes 

M ISfl) = i Jf _JL_ M<l) . (12) 

m z 7 a 3 

However, for both cases, the most conservative bound on a 3 results from maximizing 
X, minimizing c g and maximizing c^. This is accomplished by taking b 3 = —7, 
62 = —1/3, and 

aem= 12¥0 m *= 190GeV f{vM = -0-025 fj, = m H = 1 TeV . (13) 

I use a minimum value of «g in Equation (1) and a maximum value of ac in Equation 
(12) from the range 

' + JL _ 0.7 < — < -JL_ + JL + 1.4 (14) 



20a em 5a 3 a G 20a em 5a 3 

obtained from gauge coupling unification with a 3 TeV bound on the supersymmetric 
thresholds. 

The resulting bound on 0:3 is shown for the central and 1 — a values of sin 2 6w 
by wavy lines in Figure 1. To have 50 GeV < 777.1/2 < 1 TeV requires a 3 > 0.119. Note 
the differences in the approaches leading to the upper bound on 7711/2 indicated by 
solid lines and the lower bound on a 3 indicated by wavy lines. The upper bound on 
777.1/2 is obtained by maximizing the right side of Equation (1) with no bounds on the 
supersymmetric thresholds. The lower bound on a 3 is obtained by minimizing the 
right side of Equation (1) and maximizing the right side of Equation (12) assuming 
that 777.1/2, fi, vnii < 1 TeV. The ultimate reason for this difference is that the super- 
symmetric threshold corrections to sin 2 9w in the SSM are bounded from above by 
the assumption of universal soft supersymmetry breaking but are only bounded from 
below by naturalness. 
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As a simple example of modifications required in more realistic models, con- 
sider minimal supersymmetric £77(5). Although the predictions of gauge coupling 
unification depend on the unknown superheavy scales, one can extract the mass of 
the superheavy proton decay mediating triplets as a function of the low-energy cou- 
plings [To] , M. Unfortunately, this leaves no prediction for the supersymmetry 
breaking scale. However, gauge coupling constant unification, a 1 TeV naturalness 
bound, and limits on the proton lifetime can be combined to derive a lower bound 
on a 3 @. Using the value of sin 2 Q w in Equation (6) gives a bound of a 3 > 0.119 
in minimal supersymmteric £77(5). Careful study fl6| reveals that in minimal super- 
symmetric ££7(5), the supersymmetry breaking scale could be as large as 10 8 GeV 
and extensions of the minimal model further relax this bound. 

In addition to the theoretical prejudice for a naturalness bound on supersym- 
metry, this calculation offers the possibility of using coupling constant unification to 
place a bound on rai/2 in specific models. The excellent agreement of gauge coupling 
unification with the SSM, combined with the sharpening and shifting measurements 
of the low-energy couplings provide interesting speculation. On the one hand, apart 
from the highly unlikely possibility of a light gluino |17|], measurements of the strong 
coupling from deep-inelastic scattering [I8[, charmonium ||19|| , and T |20| give low 



values of a 3 which in the SSM and minimal supersymmetric ££7(5) require unnatu- 
rally high values of m\/2- On the other hand, the values of a 3 from jet shapes at LEP 
are increasing to values which ensure a low value of mi/2 in the SSM. To complicate 
matters, the value of sin 2 ^ from SLD [2I| is lower than that from LEP introducing 
more uncertainty in this critical input to gauge coupling unification calculations. 

Even though these types of conclusions can so far only be reached in very 
simple models like the SSM and minimal supersymmetric ££7(5), the dependence of 
the results on the unfolding experimental situation holds for some the excitement of 
a close horse race. The observation that low-energy limits of string models come very 
close to the simple models considered in this paper, and the ability of string theory 
to quantify gravitational corrections to the predictions of gauge coupling unification, 
offers hope of testing the physics of unified theories using precision low-energy mea- 
surements. 
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Figure Captions 



Figure 1: Solid lines indicate the upper bound on the supersymmetry breaking scale 
from gauge coupling unification in the SSM for central and 1 — a values of sin 2 9 = 
0.2317 ± 0.0005 with c=0.5. The supersymmetry breaking parameter mi/2 is below 
lOTeV for a 3 > 0.120 and is below 1 TeV for a 3 > 0.129. Wavy lines indicate 
a lower bound on the strong coupling from gauge coupling unification assuming a 
naturalness bound on the scale of superysmmetry breaking. Reasonable values of 
50 GeV < m 1/2 < 1 TeV require a 3 > 0.119. 

Figure 2: The value of 0:3, as a function of c, above which TO1/2 is bounded by 
1 TeV, 10 TeV is indicated by solid, dashed lines in the SSM for central and 1 — o 
values of sin 2 9 = 0.2317 ± 0.0005. The most conservative choice, /i, m# < raz 
corresponds to c = 0. A Monte Carlo search of the SSM parameter space finds that 
over 99% of the points considered give c > 0.5. 
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